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Abstract 

An ordered set-partition (or preferential arrangement) of n labeled elements represents 
a single "hierarchy"; these are enumerated by the ordered Bell numbers. In this note we 
determine the number of "hierarchical orderings" or "societies" , where the n elements are first 
partitioned into m < n subsets and a hierarchy is specified for each subset. We also consider 
the unlabeled case, where the ordered Bell numbers are replaced by the composition numbers. 
If there is only a single hierarchy, we show that the average rank of an element is asymptotic 
to n/(41og2) in the labeled case and to n/4 in the unlabeled case. 
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1 Introduction 

Suppose we are given a set S of n labeled elements (or "individuals"). The number of ordered 
set-partitions (or "hierarchies") on S is given by the ordered Bell number B n (see sequence 
A670 in [10] for properties and references). A hierarchical ordering or society on S is specified 
by first distributing the elements into m < n unlabeled and nonempty subsets, and forming 
an ordered set-partition in each subset. In Section 2 we will determine the number of different 
structures of this type. Section 3 discusses the analogous question in the unlabeled case, when 
the elements are indistinguishable. 

The original motivation for this work was to try to describe the structure of a "typical" 
society. Of course Sections 2 and 3 just enumerate them. However, in Section 4 we consider 
the distribution of ranks in a random selection of a single hierarchy, in both the labeled and 
unlabeled cases. 

Structures of the kind considered here were discussed in a classic combinatorics paper by 
Motzkin [6]. In his terminology, the labeled structures in Section 2 would be called "sets of 
lists of sets" , and the unlabeled structures in Section 3 "sets of lists of numbers" .* 



"The Maple package combstruct [5] makes it easy to generate such structures. For example, the labeled 



2 The labeled case 



Let H n denote the number of possible hierarchical orderings or societies, with exponential 
generating function (or e.g.f.) H(x) = Y2n>o H n x n /n\. An explicit formula for H{x) follows 
from a standard application of the exponential formula in combinatorics. 

Theorem 1 

H(x) = ex V (B(x) - 1) , (1) 
where B(x) = 1/(2 — e x ) is the e.g.f. for the ordered Bell numbers. 

Proof. This is immediate from (for example) [11, Cor. 3.4.1], or [2, Chap. 1.4, p. 46]. ■ 
The first few values H n for n = 0, 1, 2, . . . are 1, 1, 4, 23, 173, 1602, 17575, 222497, 3188806, 

50988405, . . . (this is now sequence A75729 in [10]). Table I illustrates the case n = 3. 
Several properties can be deduced from Theorem 1. 

(i) By logarithmic differentiation of (1) (cf. [11, Chap. 1, p. 22]) we obtain a recurrence 

H n = J~) B k H n _ k . (2) 

(ii) Expanding the right-hand side of (1) leads to an explicit formula: 

(mi,m 2 ,..,m„) lx J-t •> w > 

where the sum is over all (mi, m.2, • • • , m n ) such that X^j=i i m i = n > that is, m,j is the number 
of subsets of S containing j elements, for j = 1, . . . , n. 
When n = 6, for example, 

tt _ 6! I j 6! i4 ol j 6! i2 o2 , 6! o3 j 6! i3 10I , 

~~ 6!(l!) e 1 "f" 4!1!(1!) 4 (2!) 1 1 ° "T 2!2!(1!) 2 (2!) 2 1 ° 3!(2!) a ° 3!l!(l!) a (3!) 1 1 10 "T 

6! il ql ioI I 6! i q2 , 6! i2 7 cl i 6! ql 7 rl i 6! q 1 c 4 il, 

1!1!1!(1!) 1 (2!) 1 (3!) 1 1 ° io " r 2!l!(3!) 2 10 ^ 2!1!(1!) 2 (4!) 1 1 ' d ^ 1!1!(2!) 1 (4!) 1 ° ' J " t l!l!(l!) 1 (5!) 1 1 ^ 

^ 4683 1 = 17575. 

(iii) The e.g.f. satisfies the differential equation 

exp(x) 

(iv) Asymptotic behavior. The e.g.f. has an essential singularity at x = log 2. The 
asymptotics are sufficiently complicated that we sought computer assistance. Salvy's Maple 
package gdev [5], [9] is specifically designed for this purpose. The result, after simplification, 
is that 

n , e V2n/log2 

H n ~ — -j-, — -tt-, r - as n — > oo , (4) 

CV4 n 3/4 ( log 2 )n ' ^ ' 



structures in Section 2 are defined by the specification 

[H, {H = Set(Sequence(Set(Z, card >= l),mrd >= 1))}, ZafeeZed]. 

[10] now contains the sequences enumerating all such Motzkin-type structures in which the specification involves 
up to three occurrences of Set and Sequence. 
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Table I: For n = 3 elements there are 23 hierarchical orderings into at most three different 
subsets. The subsets are separated by bars, and the hierarchy within a subset is indicated by 
the vertical arrangement. 
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where 



This implies that 



C = 32 7T 2 exp(3 - — — ) log 2 = 1038.97... 
log 2 



In 

logH n ~ nlogn - nil + log log 2) + W- — - + 0(logn) 

log 2 



For comparison, B n satisfies 



Bn ~ 2(log2)«+i ' (5) 
log S n ~ ?i log n — n(l + log log 2) + 0(log n) . 
By iterating Theorem 1 we can also count hierarchies of hierarchical orderings. The e.g.f. 

is 

6XP ( 6XP G-eUr) " V " V ' (6) 

and the first few terms are 1, 1, 6, 52, 588, 8174, 134537, 2554647, 54909468, 1316675221, . . . (se- 
quence A75756 in [10]). 



3 The unlabeled case 

If the initial n elements are unlabeled, ordered set-partitions are called "compositions" of n, 
and their number is 2 n_1 [8, p. 124]. To obtain a hierarchical ordering we partition the 
elements into m < n unlabeled and nonempty subsets, and form a composition of each subset. 
Let U n denote the number of such hierarchical orderings, with ordinary generating function 
(oro.g.f.) U(x) = J2 n>0 U n x n . 



Theorem 2 

Proof. This is the unlabeled analogue of Theorem 1. If a n , n > 1, is the number of n-element 
objects with a certain property, then b n , the number of disjoint unions of such objects with a 
total of n elements, where the order of the components is unimportant, is given by 

oo oo 
n=l j=l 

(see for example [3, p. 91]). ■ 

The first few values U n are 1, 1, 3, 7, 18, 42, 104, 244, 585, 1373, . . . (sequence A34691 in [10]). 
When n = 3, for example, the H3 = 23 hierarchical orderings in Table I reduce to U3 = 7 when 
the labels are removed. 
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Properties, (i) Logarithmic differentiation of (7) leads to a recurrence: 

1 " 

U n = -^2a k U n _ k , where a k = ^d2 d ' 1 . (8) 

k=l d\k 

(ii) It is not so straightforward to find the asymptotic behavior in this case. This is to be 
expected, since the generating function for the number of partitions of n has a similar form to 
(7). Also, (7) does not belong to the family of generalized partition functions considered by 
Meinardus and discussed in [1, Chap. 7]. However, the saddle point method applies. 



Theorem 3 



2n e V2n 

U n ~ ; ; ; — CIS Tl — > CO . (9) 



Proof. We have 



00 00 on— 1 ry.mn 

\o g u{ x ) = 

^-^ — ' m 

n=l m=l 
N=l d\N 



the interchange of summations being justified since all terms are positive, so 

00 k 

l °su(x) = Ett^' ( 10 ) 

k=l 

This has poles at 2" 1 , . . . , and the radius of convergence is 1/2. The pole at 

x = 1/2 dominates, and we apply the saddle-point method as in [7, §12]. The saddle point is 
at x = r n , the solution to 

U'(r n ) 
Tn U(r n ) 

which is 



n 



1 1 / 1 

T n = --— V8n + 1 + — + 0(n 

2 on on 

Then [7, Eq. (12.9)] leads to (9). ■ 



-3/2 n 



4 The structure of a random hierarchy 

In this section we consider the case of a single (labeled or unlabeled) hierarchy. The elements 
at the bottom of the hierarchy will be said to have rank 1, those at the next level rank 2, and 
so on. The maximal rank in a hierarchy is its height. 

Suppose there are n labeled elements. Let X be a hierarchy chosen at random from the B n 
possibilities, and let x € X be a randomly chosen element. There are fo! {^} ways that X can 
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have height h, where {^} is a Stirling number of the second kind (cf. [4, Chap. 7, Problem 
44]), and indeed 



i>{;;} • 



B n 

h=0 

Given that X has height h, x € X is equally likely to have any rank from 1 to r. It follows 
that the probability that a randomly chosen x has rank r is 

The average rank is 

r=l i=l 

The numbers a n B n , n = 0, 1, . . ., have e.g.f. 

^ ™ (i+ l)! jn^i 1 (e x - l)(e x - 3) 

^~\^ 2 UJ ~ ~2 (e* - 2) 2 ' 

?1=0 1=1 V ^ 

after some simplification, using standard properties of Stirling numbers (cf. [4, Table 351]). 
This implies 

n 

o- n B n ~ n\ — : — ttt as n — > oo . (13) 

8(log2)™+ 2 v ; 

Combining this with (5) we find that the average rank is asymptotic to 



71 



4 log 2 



0.36067... 7i as n -> oo . (14) 



In the unlabeled case, if one of the 2 n 1 compositions of n is chosen at random, X (say), 
and x € X is picked at random, one can show that the probability that x has rank r is 



n2 n-l Z-^ \i 
i=r 

and the average rank is (n + 3)/4. Ranks are higher in the labeled case. 
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